Large Rapidity Gap Processes 
in Proton-Nucleus Collisions 
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The cross sections for a variety of channels of proton-nucleus interaction associated with large 
gaps in rapidity are calculated within the Glauber-Gribov theory. We found inelastic shadowing 
corrections to be dramatically enhanced for such events. We employ the light-cone dipole formalism 
which allows to calculate the inelastic corrections to all orders of the multiple interaction. Although 
Gribov corrections are known to make nuclear matter more transparent, we demonstrate that in 
some instances they lead to an opaqueness. Numerical calculations are performed for the energies 
of the HERA-B experiment, and the RHIC-LHC colliders. 
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I. INTRODUCTION 

Large rapidity gap (LRG) events in hadronic collisions 
at high energies contain important information about 
the interaction dynamics. The dominant contribution 
to these events comes from colorless gluonic exchange 
which is usually associated with diffraction or Pomeron 
exchange. The probability to have a colored exchange, 
but radiate no gluons within a LRG, attenuates exponen- 
tially as function of the rapidity gap, similar to what is 
known for secondary Reggeons. 

A colorless gluonic exchange may lead to excitation of 
the colliding proton, or the nucleus, or both. We assume 
the nuclei to be mainly a composition of colorless clus- 
ters. Therefore, a nucleus may be excited either with or 
without excitation of the clusters. The latter possibility 
is usually called quasielastic or quasidiffractive scatter- 
ing. 

The proton can also be excited via different mecha- 
nisms, with or without excitation of gluonic degrees of 
freedom. The latter possibility means that the gluons 
are not resolved by the interaction, but only the valence 
quark skeleton of the proton is excited. This channel of 
excitation greatly affects the survival probability for a 
proton propagating through the nuclear matter. Its dy- 
namics is very much model dependent. In terms of the 
Regge phenomenology, this excitation channel is related 
to the 1P1PM triple-Regge graph. 

The other possibility of proton excitation, due to 
diffractive gluon radiation, is usually associated with 
the triple-Pomeron (3JP) mechanism. It has been 
known since the 70s that the triple-Pomeron coupling 
is amazingly small, which means that diffractive gluon 
bremsstrahlung is unusually weak. This signals that the 
gluons in the proton are located within small spots. In- 
deed, large mass diffraction provides a unique access to 
the process of gluon radiation. The existence of the 1/M 2 
tail in the effective mass distribution of the diffractively 
excited hadron is a clear evidence for radiation of a vector 
particle. In order to explain the smallness of the radia- 
tion cross section one should assume that gluonic clouds 



in the proton are as small as 0.3 fm [1-4]. 

The observed smallness of of the triple-Pomeron 
diffraction in pp collisions also leads to weak gluon shad- 
owing in nuclei [1]. This is indeed confirmed by the latest 
DGLAP analysis of deep-inelastic data in the next-to- 
leading approximation [5]. 

Thus one should conclude that diffraction, in particular 
in nuclei, provides also sensitive tools for the study of the 
hadronic structure. 

Elastic scattering should be also classified as a LRG 
process. There is, however, an important difference be- 
tween elastic and inelastic diffraction, in particular, how 
much they are affected by absorptive corrections. To see 
that one can just glance at the data. While the total 
and elastic cross sections rise with energy with a rate 
corresponding to the Pomeron intercept ap(Q) « 1.1, 
the single diffraction cross section is nearly independent 
of energy [6, 7]. This may be related to the difference 
in the absorption corrections which look especially sim- 
ple in the eikonal approximation and impact parameter 
representation, 
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f sd (b)=f sd (b) e*W> 



(1) 
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Here / e ;, f s d and f e i, f s( i are the input and absorption 
corrected clastic and single diffractive amplitudes respec- 
tively. Expanding the exponentials up to the lowest order 
correction one gets a correction factor [l+£i/ e /(6)], where 
£ = 1/2 for the elastic amplitude, Eq. (1), and £ = 1 for 
the diffractive amplitude, Eq. (2). We conclude that even 
for the hadronic amplitude the absorption correction for 
the diffractive amplitude is twice as big a for the elastic 
one. Correspondingly, the result of the absorptive cor- 
rections is more pronounced in diffraction than in elastic 
scattering. In particular, the energy rise of the exponent 
in the absorption factor in (2) considerably slows down 
the energy dependence of the diffraction cross section. In 
the case of diffraction on nuclear targets absorptive cor- 
rections produce dramatic reduction of the cross sections. 
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Below we demonstrate that the probability of diffrac- 
tive excitation of the valence quark skeleton in the pro- 
ton is also quite small, only 6.5% of the elastic scattering. 
However, on the contrary to gluon radiation such a small- 
ness is due to small overlap of the initial and final wave 
functions. Therefore, this probability is sensitive to the 
proton structure as well as to the form of the dipole cross 
section. 

In the case of nuclear targets the LRG processes are es- 
pecially sensitive to the Gribov inelastic shadowing. In- 
deed, those corrections affect the nuclear transparency, 
i.e. the exponential term in (1), (2). For a heavy nucleus 
this term is tiny, therefore the inelastic corrections to the 
elastic amplitude Eq. (1) cannot be large even if the ex- 
ponential term changes considerably. At the same time, 
single diffraction and other LRG channels are directly 
affected by inelastic shadowing and may undergo consid- 
erable modification. In what follows, we perform calcu- 
lations for different LRG processes within the Glauber 
model, as well as within the light-cone dipole approach 
which incorporates the inelastic corrections in all orders 
of multiple interaction. We find that for quasielastic 
scattering inelastic shadowing make nuclear matter more 
opaque, rather than transparent. This result is at vari- 
ance with the conventional wisdom. 

There is also a practical reason to perform reliable cal- 
culations for the LRG cross sections. Many experiments 
at HERA-B, RHIC and future measurements at LHC 
have triggering systems which are set up for the central 
rapidity region and miss LRG events. For the purpose of 
normalization one has to know the part of the inelastic 
cross section covered by the trigger, i.e. the total cross 
section minus the LRG contributions. 

The paper is organized as follows. We formulate the 
light-cone (LC) dipole description of soft hadronic reac- 
tions in Sect. II. In particular, Sect. II A presents the 
phenomcnological dipole cross section fitted to data, and 
Sect. II B describes the models for the valence quark pro- 
ton wave functions used throughout the paper. 

Sect. Ill is devoted to the process of single diffractive 
excitation in pp collisions. The analysis of data per- 
formed in Sect. Ill A demonstrates a surprisingly small 
cross section of diffraction related to excitation of the va- 
lence quark skeleton of the proton. Attempting to explain 
this smallness within different models listed in Sect. IIIB 
we conclude that only a saturated shape of the dipole 
cross section maybe considered realistic. As for the mod- 
els for the valence quark distribution in the proton, the 
truth seems to lie between the two extremes, the sym- 
metric 3q configuration and the quark-diquark structure 
with vanishing diquark size. 

The triple-Pomeron part of the diffraction related to 
diffractive gluon radiation is considered in Sect. Ill C. 
This contribution is unambiguously identified in data via 
its I/M 2 tail in the effective excitation mass distribution. 
This part of diffraction is well described by the model, 
since the quark-gluon LC wave function has been fitted 
to this data previously. 



Single diffraction, as well as any off-diagonal LRG gap 
process is subject to unitarity, or absorptive corrections. 
These corrections evaluated in Sect. HID substantially 
slow down the energy dependence of diffraction. 

Unfortunately, the calculation of the double diffractive 
cross section goes beyond the employed phenomenology 
of dipole-proton cross section. It needs more refined in- 
formation about the dipole-dipole cross section. There- 
fore, we make a simple estimate of this cross section in 
Sect. HIE, based on Regge factorization. 

In Sect. IV we switch to proton- nucleus collisions and 
present the Glauber model approach to LRG processes. 
Unfortunately most of LRG channels cannot be accessed 
in this approximation which neglects all off-diagonal 
diffractive reactions. Therefore, one should include the 
Gribov inelastic shadowing corrections which are intro- 
duced in Sect. V within the eigenstate formalism. It turns 
out that at high energies the interaction eigenstates are 
the color dipolcs, thus the Gribov corrections lead to the 
color transparency effect [8] . 

In Sect. VI cross sections of coherent processes in which 
the nucleus remains intact are calculated. These include 
elastic scattering and diffractive excitation of the pro- 
jectile proton, as well as the total cross section related 
via unitarity to the elastic amplitude. We employ differ- 
ent combinations of the two models for the proton wave 
function and two models for the dipole cross section. The 
results are presented in Table I for the energy of HERA- 
B, sfs = 41.6 GeV, in Table II for RHIC, ^ = 200 GeV, 
and in Table III for LHC, y/s = 5.5 TeV 

The reactions leading to nuclear break-up, quasielastic 
scattering and diffractive excitation of bound nucleons 
are considered in Sect. VII. One can calculate the cross 
section making use of completeness. 

Besides the excitation of the proton valence quark 
skeleton, diffractive interaction can shake gluons off the 
proton. In terms of Regge phenomenology this process is 
related to the triple Pomeron part of diffraction. On the 
other hand, it causes a reduction of gluon density in the 
proton, an effect usually called gluon shadowing. These 
phenomena are considered in Sect. IX. 



II. LIGHT-CONE DIPOLE REPRESENTATION 
A. The dipole cross section 

The cross section of interaction of a qq dipole with a 
nucleon, introduced in [8] , is usually assumed to be flavor 
independent, and to depend only on the qq transverse 
separation and energy. This independence of quark flavor 
has been proven only perturbatively, but good agreement 
between data and a parameter free calculation [9] of the 
J/*!? photoproduction cross section shows that this is a 
good approximation. 

The dipole cross section calculated perturbatively [8] 
within the two-gluon model for the Pomeron [10, 11] de- 
pends only on the transverse dipole size rx and is in- 
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dependent of energy and the relative fractions of the 
momentum carried by the quark and antiquark. How- 
ever, higher order corrections bring forth a dependence 
on energy (or Bjorken x-dependence in DIS) and on the 
fractions a and 1 — a of the dipole light-cone momen- 
tum carried by the q and q. This energy dependence is 
rather mild (especially in soft processes under considera- 
tion) and may lead to sizeable effects only for large vari- 
ations of the collisions energy. Since the variation of a is 
usually not very large, its dependence can be neglected. 
This approximation may be questionable in DIS where 
the end-point region of the a distribution is enhanced by 
the increasing cross section [12]. However, this is not the 
case in soft reactions where the a-independence should 
be a good approximation. 

In what follows we will test two popular models for 
the dipole cross section: 

(i) Simple extrapolation of the small-rx be- 
havior, 



(Jqq(r T ,s) = C(s)r| ! 



(3) 



(ii) The saturated cross section rises as r^, at small 
rj., but levels off at large rf,, 



j(r T ,s) = a (s) 



1 — exp 



R 2 (s) 



(4) 



where R (s) = 0.88 fm {s /s)°- u and s = 1000 GeV 2 
[1]. The energy dependent factor cro(s) is defined as, 



0o (s) 



np 



(s) 1 + 



(5) 



where (r 2 ^ = 0.44 ± 0.01 fm 2 [13] is the mean square 
of the pion charge radius. 

This dipole cross section is normalized to reproduce 
the pion-proton total cross section, {<J qq )-K = &tot( s )- The 
saturated shape of the dipole cross section is inspired by 
the popular parametrization given in Ref. [14], which is 
fitted to the low- a; and high Q 2 data for F%(x, Q 2 ) from 
HERA. However, that should not be used for our pur- 
pose, since is unable to provide the correct energy de- 
pendence of hadronic cross sections. Namely, the pion- 
proton cross section cannot exceed 23 mb [50]. Besides, 
Bjorken x is not a proper variable for soft reactions, since 
at small Q 2 the value of x is large even at low energies. 
The s-dependent dipole cross section Eq. (4) was fitted [1] 
to data for hadronic cross sections, real photoproduction 
and low-Q 2 HERA data for the proton structure func- 
tion. The cross section (3) averaged with the pion wave 
function squared (see below) automatically reproduces 
the pion-proton cross section. 

In the case of a proton beam one needs a cross sec- 
tion for a three-quark dipole, a^ q {ri, 7=2, 7=3), where fj 
are the transverse quark separation with a condition 
7=1 + 7=2 + 7=3 = 0. In order to avoid the introduction of 
a new unknown phenomenological quantity, we express 



the three-body dipole cross section via the conventional 
dipole cross section a qq [16], 

03 9 (7=i, 7=2,7=3) = ^ 099(7-1) + cr qq (r 2 ) + (r qq (r 3 ) . (6) 

This form satisfies the limiting conditions, namely, turns 
into cr qq (r) if one of three separations is zero. Since all 
these cross sections involve nonperturbativc effects, this 
relation hardly can be proven, but should be treated as 
a plausible assumption. 



B. The light-cone wave function of the proton 

Since the dipole cross section is assumed to be indepen- 
dent of the sharing of the light-cone momentum among 
the quarks, the wave function squared of the valence Fock 
component of the proton, |$(fj, ctj | 2 should be integrated 
over fractions at. The result depends only on transverse 
separations n . The form of the nonperturbativc valence 
quark distribution is unknown, therefore for the sake of 
simplicity we assume the Gaussian form, 

|*jv(ri,f2,f 3 )| 2 



i=l 



3 = 1 



2 + r 2 p /R 2 p 
(nr p R p y 
x 5(fi +r 2 + r 3 ) 



(7) 



where n are the interquark transverse distances. The two 
scales r p and R p characterizing the mean transverse size 
of a diquark and the mean distances to the third quark. 
In what follows we consider two extreme possibilities: 

(i) The quark-diquark structure of the proton. 
There are evidences from data that the dominant config- 
uration of valence quarks in the proton contains a small 
isoscalar ud diquark of a size r p <~ 0.2 — 0.4 fm [17-19], 
i.e. r 2 <C Rp. Neglecting the diquark radius we arrive at 
a meson-type color dipole structure, 



1*^(7=1,7=2, r 3 )\ = 



irR 2 



exp - 



2 4 
R 2 



(8) 



where r=r = 7=2 = 7=3 , and R p is related to the mean charge 
radius squared of the proton as R 2 = ^§-(r 2 h ) p - 

(ii) The symmetric proton structure. Another 
extreme would be to say that the forces binding the va- 
lence quarks are of an iso-scalar nature, therefore the 



quark distribution is symmetric, i.e. r p 



R 



v in (7). 

In this case the mean interquark separation squared is 
(r 2 ) = lR 2 p = 2(r c \) p . 

Thus, we consider two possibilities for the shape of 
the dipole cross section and two choices for the proton 
wave function. In what follows we will test these models 
comparing to data. 
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III. SINGLE DIFFRACTION pp ->■ Xp 

Now we have a choice for the key ingredients of the 
dipole approach [8], which are the universal dipole cross 
section and the light-cone wave function of the incoming 
hadron. First of all, we should make sure that each model 
correctly reproduces the total pp cross section, which is 
the main entry for the calculation of nuclear effects. A 
further important and rigorous condition for a realistic 
model would be the magnitude of the single diffractive 
cross section. 



A. Experimental data 

The cross section of proton excitation, pp — ► pX, inte- 
grated over large Feynman xp > 0.95 is about 3.5 mb at 
Eiab = 920 GeV and nearly saturates at the value of 4 mb 
at higher energies [6, 7]. This cross section, however, in- 
cludes two parts: (i) excitation of the quark skeleton of 
the hadron without gluon radiation; (ii) Excitations in- 
cluding emission of gluons. The former, in terms of the 
Regge phenomenology, corresponds to the triple Regge 
graph 1P1PM and behaves like 1/M 3 at large excitation 
masses. The latter corresponds to the triple-Pomeron 
graph 31P which provides the 1/M 2 tail at large masses. 

In terms of the dipole description the lowest Fock com- 
ponent at which we are concentrating in this section cor- 
responds to the first part, P 1PM. The higher Fock com- 
ponents containing gluons will be considered in Sect. IX. 

At high energies and large xp — > 1 one can write the 
diffractive cross section in the triple- Regge form [20], 



dxp dp\ 

G 3P (0) 
(l-x F ) 



si Gpp R (0) 

7 (i-xpf 2 



-B 3B p T 



e-B&p'r + G RM p(0) e - BmBJp PT . ( 9 ) 



Here si = lGeV 2 , 



R 



2a' P ln(l — xf) 
B RRP = Rrrp - 2«k ln(l - xp) 



PPi 



(10) 



i = P,M; a' p w 0.25 GcV~ 2 and a' M w 0.9 GeV" 2 arc 
the slopes of the Pomeron and secondary Reggcon trajec- 
tories. Note that we keep in (10) only the triple-Regge 
terms which are either singular in (1 — xp), or do not 
vanish as powers of energy. 

Since the data show no substantial rise of the diffrac- 
tive cross section with energy [6, 7], which is appar- 
ently caused by strong absorptive corrections, we incor- 
porate this fact fixing the effective Pomeron intercept 
at a p (0) = 1. This also allows us to use the results 
of the comprehensive triple-Regge analysis of data per- 
formed in Ref. [20], in which a P (Q) = 1 was used. They 
found the following values for the parameters: Gsp(0) = 
G PPR (0) = 3.2 mb/ GeV 2 ; G RRP (0) = 13.2 mb/ GeV 2 ; 



H-Ip — 4.2 GeV 2 ; R PPR — 1.7 GeV 2 ; R RRP — 
0GeV~ 2 . 

We are now in a position to evaluate the diffractive 
cross section integrating the distribution Eq. (9) over p\ 
and xp, from x m i n up to x max . We fix x. m i n — 0.95 as 
is usually done to separate diffractive from nondiffrac- 
tive contributions, and x max = 1 — Mq/s with Mq = 
(mjy + m^) 2 GeV 2 . Of course such small masses do not 
satisfy the condition of the triple-Regge dynamics, but 
we appeal to the duality between s-channel resonances 
and i-channel Reggeons which is known to work well if 
one averages over the resonances. Then we get at the 
energy of HERA-B, a p sd = 3.27 which agrees well with 
the value <j p sd " — 3.5 mb suggested by data. We expect 
a sd ~ ^ m b a t R-HIC and LHC energies. 

As was already mentioned, in this section we are only 
interested in the part of diffraction related to excitation 
of the valence quark skeleton without gluon radiation. 
This part corresponds to the first term, P1P1R, in Eq. (9). 
Integrating it over p\ and xp we get a^(lPlPlR) = 
1.13 mb, i.e. about the third of the total single diffractive 
cross section. 

For further model tests we need to know the forward 
diffraction which turns out to be quite small with respect 
to the forward elastic cross section, 



Rsd 



d<r sd /dp, 



da e i/dp 2 j 



PT—O 



5.5 mb/ GeV 2 
84.5 mb/ GeV 2 



= 0.065 



(11) 

at the energy of HERA-B. 

It is worth reminding that we rely on the hypothesis 
of duality and should not expect this result to be very 
accurate. It just shows the scale of the effect, namely, 
that diffractive excitation of the projectile valence quark 
system is an order of magnitude weaker than the elastic 
channel. 

In hadronic representation, in particular within the one 
pion exchange model, one can explain rather well such a 
weak diffractive excitation of the valence Fock component 
of the proton [20]. It is a challenge, however, to reproduce 
this result within the QCD based dipole model. 



B. Models 

Since we have models for the dipole cross section and 
the proton wave function, it is therefore possible to calcu- 
late the total and forward single diffractive cross sections 

[1, 8], 



°t5t = Kg (»•<)> 



(12) 



dM 2 d(J sd(pP^pX) 



dM\dp 2 T 



p T =0 



(13) 
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Here the averaging weight is the proton wave function 
squared, 

(•••>= ffldWnr^ 2 (.••)• (14) 
J i=i 

We are going to test the models mentioned before regard- 
ing their ability to explain the strong suppression of the 
diffractive relative to the elastic channels, R sc i <~ 0.1 in 
Eq. (11). 



separation grossly overestimate the diffractive cross sec- 
tion, since the main contribution comes from large dis- 
tances. Therefore, both models I and II are not realistic 
and should not be used for comparison with data. 

Actually, Eq. (13) shows that diffraction is given by 
the dispersion of the dipole cross section dependence on 
the transverse qq separation. Apparently, (Jq q (r T ) oc re- 
used so far varies with rr too steeply and the dispersion 
is too large. One needs a flatter dependence in order to 
suppress diffraction. 



1. Model I: quark-diquark proton and <Jqq <x r\ 



3. Model III: quark-diquark proton and saturated cross 
section 



The simple quadratic tt dependence of the dipole is 
quite popular in the literature and is frequently used for 
calculations of nuclear effects (e.g. see [16, 21, 22]). Since 
the energy dependent factor is unknown, it is usually fit- 
ted to reproduce the total cross section. In this case 
one has to explain the relative values of total cross sec- 
tions for different hadronic species, in particular, the ra- 
tio afo t / '<jJo t w 1.6. Indeed, the ratio of charged radiuses 
squared agrees with this value within the diversity of the 
results of different measurements [23]. 

However, the cross section of single diffraction turns 
out to be dramatically overestimated. Indeed, the diffrac- 
tion to elastic ratio defined in (11) reads, 



Rsd 



(4) 
(4> 2 



-1 = 1. 



(15) 



This is more than one order of magnitude larger than the 
experimental value Eq. 11. 



2. Model II: symmetric proton and (Tqq oc r% 

In this case we again treat the unknown factor in the 
dipole cross section as a free parameter and adjust it to 
the total cross section, 

(vz q m = vZ ; (is) 

Then we can calculate the mean value of the dipole cross 
section squared, 



3 



(17) 



and eventually the diffractive cross section according to 
(13). Then we arrive at the ratio, 



1 



Rsd — 



which is closer to the experimental value Eq. (11), but 
still quite overestimates the data. 

It is not a surprise that both above models based on 
the dipole cross section Eq. (3) steeply rising with quark 



There is no freedom in this case, the dipole cross sec- 
tion Eq. (4) if fixed to reproduce the pion-proton cross 
section. The calculated value of proton-proton total cross 
section is rather close to the experimental value, but 
somewhat smaller. This is not a surprise, since this 
model employs the approximation of a point-like diquark, 
which apparently leads to an underestimation of the pp 
cross section. It is risky, however, to rely on such an ap- 
proximation for calculating nuclear effects. Even a small 
deviation of o v t ^t from the experimental value will cause 
nuclear effects which may be misinterpreted as inelastic 
shadowing. To make sure that this cross section is ex- 
actly reproduced, we redefine the function cr (s) in (5), 



(19) 



where the parameter S — 8(r 2 h ) p /3i? 2 (s), and we use the 
mean proton charge radius squared (r^ h } p = 0.8 fm 2 [23]. 

Using the proton wave function Eq. (8) and the satu- 
rated cross section Eq. (4) we get the single diffractive 
cross section in the form, 



dp% 



PT—O 



S 2 

(1 + £) 2 (1 + 26) 16tt 



Dividing by the elastic cross section we get, 



1 



1 + 26 



(20) 



(21) 



This ratio is rather small at the energy of HERA-B, 
R s d = 0.13, which is compatible with the data. The 
fraction of single diffraction decreases with energy down 
to R sd = 0.06 at the energy of RHIC {y/s = 200 GeV), 
and R sd = 0.01 at the energy of LHC (y/s = 5.5 TeV). 

One may wonder why the saturated cross section leads 
to such a weak diffraction and why decreases with en- 

(18) 

ergy? This is easy to interpret. Indeed, if the cross 
section were completely flat, i.e. cr^r^) — const, no 
diffraction would be possible because of orthogonality of 
the initial and final valence quark wave functions. Only 
the drop of <7g g (rr) at tt < Ro{s) makes diffraction pos- 
sible. However, Rq(s) decreases with energy, therefore 
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the shape of the dipole cross section is getting flatter 
and diffraction vanishes. As far as diffraction gets its 
main contribution from small rr < Ro, note that it is 
less probable to find three quarks with small separations, 
than a two-quark system. Therefore, one should expect 
less diffraction for a symmetric three-quark wave function 
of the proton as is demonstrated in the next section. 



energies much higher than those where data for <r^(s) 
are available. Then the pion to proton ratio of the cross 
sections reads, 



^toti s ) 



JPP 
'tot 



l Rlis) 



3 ( r ch)p 



(*) 1 + 



3 gjUf) 

8 kj: 



(26) 



4. Model IV: symmetric proton and saturated cross section This ratio slowl y rises with energy. At the energy of 



HERA-B R. 



Using the wave function Eq. (7) with r p = R p and the 
cross section (4) we get the following forward elastic cross 
section and diffraction to clastic ratio, 



dp, 



7 



Pt=0 
R.id 



(l + § 7 ) 2 167T ' 

(l + |7) 2 -i 



(l + I 7 )(l + f 7 )(l+ 7 ) 

The parameter 7 is related to previously introduced S 



(22) 
(23) 



7 = 3 



V chip 

Rl{s) 8 



9 



(24) 



Then, for the energy of HERA-B we arrive at a very small 
fraction of the single diffractive cross section R sc i = 0.07 
which is in excellent agreement with the experimental 
result Eq. (11). We expect a substantial reduction of 
this fraction at higher energies, R s( i w 0.03 at the energy 
of RHIC, and R sd w 0.0045 at the energy of LHC. 

Note that the the parameter 7 in Eqs. (22)-(23) can be 
defined differently, 



'tot 



(<73g) = CO (s) 



7 



l + §7 



(25) 



where ao(s) and CT3 g (fj) are defined in (5) and (6) respec- 
tively. 

Both definitions (24) and (25) nicely agree at the en- 
ergies of fixed target experiments. More problematic is 
to apply Eq. (25) at high energies of colliders, RHIC and 
LHC. No data for a^ t are available at energies above 
t/s = 30 GeV. The usual extrapolation with a universal 
energy dependence as in pp collisions is just an educated 
guess not supported by any dynamic theory. Moreover, 
one should expect a steeper rise of a^ t than af^ t . Indeed, 
the rising part of the cross section related to gluon ra- 
diation is nearly independent of hadronic size, while the 
constant part of the cross section related to dipole-dipolc 
collision followed by no gluon radiation apparently de- 
pends on the dipole sizes and is smaller for np, than for 
pp. Thus, pion-proton cross section rises steeper with en- 
ergy, in accordance with the general trend of steeper en- 
ergy dependence for smaller dipolcs discovered at HERA. 

One can combine Eqs. (24) and (25) in order to find 
the unknown energy dependence of <r£^(s). This is a 
more reliable procedure, since the dipole cross section, 
in particular the parameter Rq(s), is fitted to data at 



ir/p 



0.6 in excellent agreement with data, 



but it becomes 10% larger, R w / P — 0.66, at the energy 
of LHC. Eventually, at very high energy when unitarity 
will be saturated, all the cross sections reach the universal 
Froissart limit corresponding to an expanded black disk. 
Correspondingly, R v / p (s) — > 1. 



C. Diffractive gluon radiation 

Besides excitation of the valence quark skeleton of the 
proton, a valence quark itself can be excited followed by 
gluon radiation. In terms of Regge phenomenology this 
process corresponds to the triplc-Pomcron contribution 
to the diffraction cross section. This is the second term 
in Eq. (9). 

It has been known since the 70s that the triple- 
Pomeron coupling is quite small. To appreciate this 
statement one can express diffraction in terms of the 
Pomeron-proton total cross section which should be ex- 
pected to be twice as large as a meson-proton one. In- 
deed, the Pomeron is a gluonic system, therefore one 
should have an extra Casimir factor 9/4 compared to a 
quark dipole. Thus, one may expect crf^f s=s 50 mb. This 
expectation is in dramatic contradiction with data [6] 
which show only <rf^f = 2 mb at large excitation masses. 
Apparently, it is much more difficult to shake gluons off 
valence quarks, compared to pQCD expectations. The 
way out of this puzzle is to suggest that gluons in the 
proton are located within small spots which are hardly 
resolved by soft interactions. The mean transverse size 
of these spots was fitted to single diffraction data with 
large effective masses and found to be (r-r) ~ ro = 0.3 fm. 
Such gluons have much more intensive Fermi motion than 
massless perturbative ones, and they are less sensitive to 
an external kick, i.e. gluon radiation is suppressed. 

Such a picture is quite successful in explaining the en- 
ergy dependence of the total hadronic cross sections and 
elastic hadronic slopes [2]. It also correctly predicts the 
reduced value of a'p, the slope of the Pomeron trajec- 
tory for the process of elastic photoproduction of J/^ 
[3]. Similar statements about gluonic structure of the 
proton has been done in [4] recently. 

The light-cone wave function of the quark-gluon Fock 
component of a quark was calculated in [1] within a model 
describing the nonperturbative interaction of gluons via 
a phenomcnological light-cone potential taken in an os- 



7 



cillatory form. The result reads, 

= ~—\hi 12— 



exp 



2r 2 



^ • (27) 



Here we assume that the gluon, which is a vector parti- 
cle and possesses a l/i distribution, is carrying a neg- 
ligible fraction (x) C 1 of the quark momentum. Of 
course the concrete shape of the light-cone potential is 
not crucial. What is only important is the smallness of 
the mean quark-gluon separation. Notice that Eq. (27) 
can be viewed as a source of small gluonic spots in the 
proton [3]. 

Since the light-cone quark-gluon distribution function 
is known, one can calculate the cross section of diffractive 
gluon radiation by a high-energy quark interaction with 
a nucleon [1], 



da(qN -» qGN) 



dx f dp\ 



d 2 r T 



167r(l — xp) 
a z (r T ,s). (28) 



Pt=0 
2 

2 



Here the cross section a(rT,s) is not just a cross sec- 
tion of interaction of a qG dipole. This dipole is not 
even colorless. As usual, diffractive excitation is possible 
due to a difference between elastic amplitudes for differ- 
ent Fock states, in this case the bare quark \q) and the 
\qG) pair. Since they have the same color, the difference 
emerges from the color-dipolc moment of the q ~~ G pair. 
It is shown in [1] (see in particular Appendix A. 2) that 

v{r T ,s) = §£7g g (r T ,s). 

The next step is to integrate over px the cross section of 
diffractive gluon radiation provided that the forward one, 
Eq. (28), is known. In terms of the Regge phenomenology 
diffractive radiation corresponds to the triple-Pomeron 
term in the cross section of single diffraction. Data agree 
with a Gaussian pr-dependence of the triple-Pomeron 
term with the slope [20] , 



B p 3P (x F ) = B° 3P + 2a' P ln ( T ^—) 

\ 1 — XF J 



(29) 



where B° p 



4.2 GcV~ 2 , and a' p = 0.25 GcV~ 2 . 
Now, we are in a position to evaluate the effective 
triple-Pomeron part of the single diffraction cross section 
for pp collisions employing the the wave function Eq. (27) 
and the saturated cross section Eq. (4), 



da(pp — > pX) 



x In 



dxF dp T 
, e 2 ( S ) 



l + 2e(s) 



81a s q 2 (s) 
3JP (16tt) 2 (1 - x F ) 

exp[-p T B p 3P (x F )] , (30) 



where e(s) = r 2 , / Rq (s) . In the energy range under discus- 
sion e(s) is rather small, then the single diffractive cross 
section Eq. (30) is proportional to r$. This is why this 
process is quite sensitive to the value of t-q and provides 



an efficient way to determine the size of gluonic spots in 
the proton, ro w 0.3 fm [1]. 

Notice that the interference between the diffractive am- 
plitudes of gluon radiation by different quarks in the pro- 
ton should not be appreciable, since ro is small compared 
to the proton radius. Explicit calculations performed in 
[1] confirm this. 



D. Unitarity corrections 

Any large rapidity gap process is subject to unitarity or 
absorptive corrections which may be substantial. Indeed, 
initial/final state interactions tend to fill the rapidity gap 
by producing particles, and one may treat such correc- 
tions as a survival probability of the rapidity gap. Such 
corrections become especially large and may completely 
terminate the gap in the vicinity of the unitarity limit, 
which is also called black disk regime. 

Since the phenomenological dipole cross section <Jq q (s) 
is fitted to data, we assume it to include by default all 
the unitarity corrections. The same is true for the off- 
diagonal amplitude of diffractive excitation of the valence 
quark system without gluon radiation, since it is a lin- 
ear combination of elastic dipole amplitudes. Thus, our 
calculations for the JPIPM term in (9) do not need any 
unitarity corrections. 

The situation is different for the WIPP term, or gluon 
diffractive radiation. Although the amplitude of diffrac- 
tive excitation of a quark, qN — ► qGN, does include all 
the absorptive corrections contained in the phenomeno- 
logical dipole cross section, the presence of other projec- 
tile valence quarks, the spectators, should not be ignored. 
Indeed, any inelastic interactions of the large-size three 
quark dipole in the proton, will cause particle produc- 
tion which will fill the rapidity gap. Thus, one may ex- 
pect large absorptive corrections to the cross section of 
diffractive gluon radiation. 

Data for elastic pp scattering show that the partial 
amplitude f^(b,s) hardly varies with energy at small 
impact parameters b — > 0, while rises as function of en- 
ergy at large b [2, 24, 25]. This is usually interpreted 
as a manifestation of saturation of the unitarity limit, 
Im/J < 1. Indeed, this condition imposes a tight re- 
striction at small b, where Im/J 1, leaving almost no 
room for further rise. Correspondingly, the amplitude of 
any off-diagonal process including single triple-Pomeron 
diffraction acquires a suppression factor 

f!S(b, s) => f%(b, s) [1 - Im/J(6, a)] , (31) 

due to unitarity or absorptive corrections. This factor re- 
lated to the survival probability of LRG is known to de- 
crease with energy [26] . Interply of the rising and falling 
energy dependences of the two factors in (31) may ex- 
plain the observed flat behavior of the single diffractive 
cross section [6, 7]. 

Since Im f p f{b, s) is known directly from data, it would 
be straightforward just to fit the data with any proper 
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parameterization and use the result in Eq. (31). Alter- 
natively, one can use any model which provides a good 
description for lmf pp (b, s). It is demonstrated in [2] that 
even the simple model treating the Pomeron as a Regge 
pole with no unitarity corrections describes reasonably 
well not only the total hadronic cross sections, but even 
the data for f p [(b, s). In this model, 



Im/J(M 



(s) 



4nB PP (s) 



exp 



2B p J(s) 



(32) 



Here and for further applications we use the parametriza- 
tionfrom [27], a pp t {s) = 18.76 mb x (s/M 2 ) - 093 + a pp (s), 
where M = 1 GeV, and the Reggeon part of the cross 
section c^ p (s) is small at high energies and can be found 
in [27]. The elastic slope B pp {s) = B% + 2a' P ln(s/M 2 ) 
with B a el = 8.9 GeV" 2 and a' P = 0.25 GcV~ 2 . Note that 
the elastic amplitude at small impact parameters, i.e. 
the pre-exponential factor in (32), hardly changes with 
energy imitating saturation of unitarity. This fact has 
been known back in the 70s as a geometrical scaling. It 
is demonstrated in [2] (see Fig. 9) that not only at b = 0, 
but in the whole range of impact parameters the model 
Eq. (32) reasonably well describes the energy dependence 
of the partial amplitude f p [(b, s). 

Using Fourier transformed Eq. (32) we arrive at the fol- 
lowing cross section for single diffraction integrated over 
momentum transfer, 



da(pp — ► pX) 



d,XF 



81 a s cr 2 (s) 



3iP 



x In 



1 + 



* 2 (*) 
l + 2e(s) 



1 - 



x F ) B P p (x F ) 
nB p p (x F ) + 2B pp (s) 



(16tt) 2 (1 
1 



(47T) 2 B%(s) [B pp p (x F )+B pp { S ) 



At the scale corresponding to the mean transverse 
momentum of gluons, a s (l/r 2 ) w 0.4 [2], (33) agrees 
with data reasonably well. Nonetheless, we think that 
one should perform more elaborated calculations when 
Im/ e ;(s) is close to the unitarity limit. Indeed, the ab- 
sorptive correction factor in (31) is so small in this case 
that any little variations of Im f e i (s) may lad to dramatic 
effects. Such a fine tuning goes beyond the scope of the 
paper, but we plan to work more on this elsewhere. 



Besides, neither perturbative QCD calculations [28], nor 
phenomcnological dipole cross sections fitted to DIS data 
[14] confirm Regge factorization. The very Q 2 depen- 
dence of the Pomeron intercept observed at HERA is a 
direct evidence of lack of factorization. 

Unfortunately, the phenomenological dipole cross sec- 
tion (jq q which is averaged over the size of the target pro- 
ton, is not sufficient for calculating double diffraction. 
Therefore, we will employ the approximate relation (34) 
in what follows, hoping that the corrections are not much 
larger than in other known cases. 

For the integrated double diffractive cross section re- 
lation (34) can be rewritten as, 



pp _ 
"dd — 



„pp 



(B P S 



p pp R pp 

a el a dd 



(35) 



At the energy of HERA-B the elastic slope B p f = 
12.6 GeV" 2 ; a pp = (a pp ot f / (l^B pp ) = 6.7mb; o% = 
3.5 mb. The slope of the IP 1PM part of single diffrac- 
tion relevant for Eq. (35), varies with Feynman x F , 
B Z = B o - 2a 'p M 1 - x F ), where a'p = 0.25 GeV~ 2 , 
B « 2 GeV" 2 [20]. The mass distribution of diffractive 
excitation of the valence quark skeleton of the proton 
strongly peaks in the resonance region at M w 1.5 GeV, 
which corresponds to 1 — x F = M 2 /s. Then B vp d = 
5.3 GeV -2 . As for the double diffractive slope, it con- 
tains only the Pomeron contribution, B pp d — —2a' p ln(l— 
x F ) = 3.3GcV~ 2 . Eventually, we arrive at an estimate 
for the double diffractive cross section at the energy of 
HERA-B, a pp = 1.25 mb. The same Eq. (35) leads to 
the double diffraction cross sections a pp . 



1.18 mb and 

a dd = 0.47 mb at the energies of heavy ions at RHIC, 
(33) sfs = 0.2 TeV, and at LHC, ^JTs = 5.5 TeV, respectively. 



IV. pA COLLISIONS: GLAUBER MODEL 

The pA clastic amplitude at impact parameter b has 
the eikonal form [29], 



A 

T pA (b; {*,, z,}) = 1 [] I 1 TPN (b- Sk) 



(36) 



fc=i 



E. Double diffraction pp — ► XY 

If the Pomeron were a true Regge pole, one would ex- 
pect Regge factorization relating the forward single and 
double diffractive cross sections, 



V dt J t=0 \ dt J t=0 \ dt J t=Q ■ 



(34) 



However, even within the Regge model this relation is 
broken due to presence of Regge cuts. The usual accu- 
racy of relations based on Regge factorization is 10—20%. 



where {sj,Zj} denote the coordinates of an i-th target 
nucleon; iT pN is the elastic scattering amplitude on a 
nuclcon normalized as, 

a P N t =2jdHReT pN (b). (37) 

In the approximation of single particle nuclear 
density[51], 



A 

pA(h,z 1 )= / n^l^MIr,-})! 2 , 

J i=2 



(38) 
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the matrix element between the nuclear ground states 
reads, 



1 - 



1 - j J d 2 sT pN (s) J dzp A (b-s,z) 



(39) 



Correspondingly, the total pA cross section has the form, 



<r^ = 2Re / d z b{l 



exp 



aC t (l-ia pp )TH(b) 



(40) 



where a pp is the ratio of the real to imaginary parts of 
the forward pp elastic amplitude; 

T » = 4v f d " s Rer pN (s)T A (b- s) ; (41) 

/oo 
dzp A (b,z) , 
-oo 



and 



(42) 



is the nuclear thickness function. We use the Gaussian 
form of T pN (s), 



ReT pN (s) 



P N 



4nB 



P N 



exp 



2B 



P N 



(43) 



Note that the accuracy of the optical (exponential) ap- 
proximation in (40) is quite good, <~ 10~ 3 for heavy nu- 
clei, but for numerical calculations we relay on the exact 
Glauber expressions throughout the paper. In what fol- 
lows we neglect the real part of the elastic amplitude 
which gives quite a small correction, ~ p 2 p /A 2 ^ 3 , and 
can be easily implemented. 

We performed numerical calculations at Ei a f, = 
920 GeV, having in mind the needs of the experiment 
HERA-B at DESY. At this energy we used a pp ot = 
41.2 mb, B PP = 12.6 GcV~ 2 and a pp = (o pp ot ) 2 /16ttB pp = 
6.7 mb. The Wood-Saxon form was used for the nuclear 
density with parameters fixed by data on electron-nucleus 
clastic scattering [30], except carbon whose density was 
taken in an oscillatory form [31], 



l>v[n : 1 t) 



1+ ar 2 } e-° r ~ 



(44) 



The parameter a was fitted to data for electron-carbon 
scattering, a = 0.0143 GeV 2 , and we assumed identical 
distributions for protons and neutrons. The results are 
depicted in Table I. 

We also performed calculations at the energies of RHIC 
and LHC, y/s — 5.5 TcV using the input cross sections 
listed at the end of Sect. HIE. The results are presented 
in Tables II and III respectively. 



A. Elastic cross section 

The simplest process with a large rapidity gap (LRG) 
is elastic scattering. It worth noting, however, that 
this channel is enhanced by absorptive corrections, while 
other LRG processes considered below are suppressed by 
these corrections. 

The elastic cross section according to (39) reads, 



pA _ 



= d l b 



1 — exp 



P N 



(45) 



Numerical results are shown in Table I. 



B. Diffractive excitation of the beam and target 

Diffractive excitation of the beam and/or target is an- 
other example of a LRG process. Unfortunately, the 
Glauber model is a single channel approximation and 
cannot treat properly diffractive excitation of the beam. 
Nevertheless, the cross section of diffractive excitation of 
the nucleus can be calculated, provided that the elastic 
and single diffractive cross sections for NN collisions are 
known. 

If the nucleus is excited without new particle produc- 
tion, i.e. it breaks up to nucleons and nuclear fragments, 
such a process, pA — ► pA* , is called quasielastic scatter- 
ing. Summing over final states of the nucleus \F) and 
applying the condition of completeness, one gets, 

771 J 



x (F \T pA {b)\0) - |(0|L pA (6)|0)| 



d 2 b 



\r pA (b)\ 



|(o|r pA (&)|o) 



.(46) 



Here the cross section of elastic pA — > pA scattering, 
Eq. (45), is subtracted. 

The fist term in this expression contains the quadratic 



term Jd 2 sT A {b - s) [T pN (s)] = T%{b)o el 
quasielastic cross section gets the form, 



P N 



Then the 



o%(pA^pA*) 

Jd 2 b {exp [-o£TZ(b) 



^ J 



d 2 bT A (b) exp 



—a: 



— exp 

pN rph 



pN T h 



mb) 



Hi 

(47) 



Another possibility to excite the nucleus is to excite a 
bound nucleon, pA — ► pY. To specify the terminology, 
following [32] we call this channel target single diffraction 
(tsd). Since the nucleus breaks up anyway, and the debris 
of the bound nucleon stay in the target fragmentation 
region, they cannot fill the rapidity gap, therefore their 
final state interaction do not affect the LRG cross section. 
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Then it must have the same form as the quasielastic cross 
section Eq. (47), except the normalization factor, 



■pY) = 



T pp 

' sd pA 
JPP qel 



(48) 



We fix the single diffractive pp cross section at a 99 , — 
3.5 mb for £/ & = 920 GeV in accordance with [6]. For the 
energies of RHIC and LHC we extrapolate the Tevatron 
data [6] assuming no energy independence, a P d = 4 mb. 

The numerical results for the target dissociation cross 
sections arc shown in Tables II and III. 



COLOR TRANSPARENCY 



The light-cone dipole representation was proposed in 
[8] as an effective tool for calculation of hadronic cross 
sections and nuclear shadowing, relying on the observa- 
tions that color dipolcs arc the eigenstates of hadronic 
interactions at high energies, and the cigenstate method 
[33] is a powerful tool for summing up all Gribov inelastic 
corrections. 

The key ingredient of this approach, the cross section 
of the dipole-nucleon interaction, <7q q (rj>), is an univer- 
sal and flavor independent function which depends only 
on transverse separation tt and energy. Applications of 
the dipole formalism to nuclei are especially simple, if 
the energy is sufficiently high to freeze the fluctuations 
of the dipole size during its propagation through the nu- 
cleus. Otherwise one should rely on the light-cone Green 
function technique [34-36] , which takes care of these fluc- 
tuations. 

Due to color screening colorless point-like dipoles can- 
not interact with an external color field. Since the under- 
lying theory is non-abelian, the interaction cross section 
for such dipoles vanishes at tt — > as <Jq q {rT) cx t\ [8], 
the phenomenon called color transparency [52]. At high 
energies nuclei are transparent for small-size fluctuations 
of the incoming hadron, therefore the exponential atten- 
uation suggested by the eikonal Glauber formula cannot 
be correct and should underestimate the nuclear trans- 
parency, i.e. overestimate the total hadron-nucleus cross 
section. 

Thus the Glauber approach is subject to modifications 
called Gribov's corrections. Originally those corrections 
were proposed in hadronic representation where they look 
like intermediate diffractive excitations [37]. The low- 
est order correction is expressed via the single diffraction 
cross section (hN — > XN) measured experimentally 



[37, 38], 



Aat A t = -4tt / d 2 6exp 



\°t N t T A (b) 



I 



dM 2 d(Tsd 



oc 



Pt- 



dM 2 dp\ 
x / dz 2 p A (b,z 1 )e lqL(z2 - Zl) 



oo 

/ dz 1 p A (b,z 1 ) 
o J 



where 



M 2 



QL = 



2E h 



(49) 



(50) 



Therefore, one might think that this is a model- 
independent calculation. However, the cross section of 
interaction of the intermediate inelastic state X is un- 
known and is assumed to be equal to . 

Formally, all the obscrvables calculated cither in color- 
dipole, or hadronic representations must be identical. 
Nevertheless, although the correction Eq. (49) is nega- 
tive, i.e. it makes the nuclear matter more transparent, 
it cannot reproduce the effect of color transparency [39] 
which results from compensation of many diagonal (pos- 
itive) and off-diagonal (negative) amplitudes. Indeed, for 
heavy nuclei this correction vanishes exponentially with 
nuclear thickness T A (b). 

The dipole representation allows to sum up the in- 
elastic corrections to all orders [8]. For a dipole of a 
fixed size tt the eikonal form is exact, since the dipole 
is the eigenstate of interaction. Therefore nuclear trans- 
parency, which is the no-interaction probability of prop- 
agation of a qq dipole with fixed separation tt through 
nuclear medium, reads, 



Tr(r T ) = e - <T « (rT)TA . 



(51) 



In averaging over dipole sizes important contribution 
comes only from small <Jq q {rT) < 1/T A . Therefore, 
for a sufficiently long path in the nuclear medium only 
very small values of tt contribute, and any model for 
the dipole cross section must have the same behavior 
Vqq{r T ) oc r%. Then [8], 



Tr = 



e -<r, Q (r T )T A ^ ^ 



1 

T~ A 



(52) 



This result should be compared with the exponential at- 
tenuation in the Glauber model, Tr = exp(— ofj^ T^). 
Such a difference cannot result from the lowest order in- 
elastic correction Eq. (49) which has the same exponen- 
tial dependence on T A . This is a full color transparency 
effect which must include all the higher order inelastic 
corrections. It is characterized by the nuclear satura- 
tion scale Q 2 oc T A and can be treated perturbatively for 
sufficiently heavy nuclei or at very high energies. In re- 
ality, such large nuclear thicknesses are beyond the reach 
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of existing nuclei. For this reason, one should not rely 
on the limiting behavior Eq. (52), but employ realistic 
models for the dipole cross section. In what follows we 
demonstrate how important is this fact. 



VI. pA COLLISIONS: THE TOTAL, ELASTIC 
AND SINGLE DIFFRACTIVE CROSS SECTIONS 

A. Excitation of the valence quark skeleton of the 
proton 

Once models for the proton wave function and the 
dipole cross section are chosen one is in a position to 
perform calculations for Gribov's corrections to different 
nuclear reactions, in particular LRG processes. The total 
and elastic cross sections read, 



of* = 2 J d 2 b [l-^e-W'-OTim 

r 3 
2 J d 2 b i-Jf[ePn |*jv(r,)| 2 



x e 2 



pA _ 



\<y 3q (r k )T A {b) 



'el 



d 2 b 



1 



■W(ri)Tji(6) 



(53) 



(54) 



Single diffractive excitation of the projectile proton, 
pA — > XA cannot be treated properly within the single- 
channel Glauber model approximation which assumes 
that the projectile hadron is an eigenstate of the interac- 
tion. To generalize the model one should introduce off- 
diagonal diffractive amplitudes, but then one faces the 
same problem as in the case of higher order Gribov's 
corrections: lack of experimental information for those 
amplitudes. And the remedy is the same, to switch to 
the dipole representation. 

The cross section of single diffraction on a nucleus 
related to excitation of the valence quark skeleton, 
[&sd]ppm, is given by Eq. (13), but with the replace- 
ment, Gqq => af q = 2 fd 2 b {1 - exp[-±aq q T A (b)}. 



dM 2 x 



= d z b 



da sd (pA -> XA) 



dM\ 

-<T3«,(ri)T*(6) 



PPP. 



e 2 



I 



(55) 



Both terms in this expression are vanishingly small for 
heavy nuclei, except at the very periphery. Therefore, 
the cross section of single diffraction is expected to rise 
as A 1 / 3 , but the coefficient should be sensitive to the 
inelastic shadowing corrections. 



Although comparison with single diffraction performed 
in previous sections for four different variants of the 
dipole model clearly demonstrated that only the last two, 
which employ the saturated dipole cross section, may be 
realistic, we will try all four cases for nuclear reactions 
to get an idea of the theoretical uncertainties. 

The inelastic corrections to of^ have been well de- 
tected experimentally [40, 41] and found to be rather 
small, about 10%. However, one should not think that 
inelastic shadowing is always a weak effect. In fact, for 
heavy nuclei it affects dramatically the exponential term 
in (40), but the term itself is very small compared to 
one, since the unitarity limit is almost saturated. How- 
ever, LRG channels (except elastic scattering), e.g. the 
single diffraction Eq. (55) , may be very sensitive to these 
corrections, since their cross section is proportional to 
nuclear transparency. 

It was demonstrated in Sect. V that for hadrons heavy 
nuclei are much more transparent than the Glauber 
model predicts, due to color transparency and the pres- 
ence of small-size dipoles in the hadronic wave func- 
tion. As a result, the exponential attenuation switches 
to an inverse linear dependence on Ta, Eq. (52). This 
is, however, an asymptotic behavior valid in the limit of 
Tjid 3 ) ~~ * 00 7 since for real nuclei the result depends on 
the actual modeling of ^Ar(^) and (Tq q (rT)- 
Model I 

Here we employ the quark-diquark model Eq. (8) for 
the proton wave function and the dipole cross section 
Cqqirr) oc r\. Then the averaged exponential terms in 
(53)-(55) read, 



exp 



1 



1 



l + |a. 



NN 
tot 



T h A {b) 



(56) 



The results of computing the total and elastic pA cross 
sections are compared with the Glauber model in Ta- 
ble I. The effect of inelastic shadowing is rather large, 
in fact the total cross section is reduced by about 20%. 
Although no data are available at this energy, extrapo- 
lation from lower energies [40] is hardly compatible with 
such a correction. Apparently this is closely related to 
the found overestimation of diffraction by this model of 
the dipole cross section. 
Model II 

Since the probability to be found in a point like configu- 
ration is less for three- than two-body system, one should 
expect more opaque nuclei in this variant. Indeed, 



/ 



exp 

3 



~^ q {n)T h A {b) 



Hd\ |vMn,r 2 ,r 3 )| 2 e -W(ft,W)Tj(i) 

i 

l 



(57) 



In this case the nuclear transparency is quadratic, rather 
than linear function of the inverse nuclear thickness [16]. 
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The results of the calculation, depicted in Table I, show 
that the inelastic correction is about 10% of the total 
cross section for heavy nuclei. This looks much better 
than for the previous model, and is compatible with what 
may be expected as an extrapolation of data at lower en- 
ergies. However, it is too early to jump to conclusions: 
the triple-Pomeron part of shadowing has not been in- 
cluded yet. 

Naturally, the elastic cross section follows the total 
cross section and is reduced by the inelastic corrections as 
well. As we mentioned, it may be at variance with other 
LRG channels which one may expect to be enhanced, if 
nuclei are more transparent. However, further calcula- 
tions show that the situation is more complicated. 
Model III 

The steep rise with r<r of the dipole cross section, 



Ti 



used above is justified only for small, but 



not for large tt- This is why it overestimates diffraction. 
More reliable calculations can be done using a more re- 
alistic phenomenological cross section which levels off at 
large rx, as described in Sect. II A. 

We can perform analytic calculations with the satu- 
rated cross section Eq. (4) expanding the Glauber expo- 
nentials. Then the mean value of the exponential terms 
in Eqs. (53)- (55) for the total cross section gets the form, 



exp 



\a qq (r T )T^b) 



ill 



= e 



-W«)t» \- Ms)r»]" 

2™ (l + n6)n\ 



E 

71=0 



(58) 



where cto(s) and 5 are defined in (19). 

As could be expected, the numerical results shown in 
Table I demonstrate a weaker effect of inelastic shadow- 
ing compared to the previous variants. 
Model IV 

This case involves the symmetric 3q-wave function of the 
proton and the saturated cross section Eq. (4). Expand- 
ing again the Glauber exponential and performing part 
of the integrations analytically, we arrive at the following 
result, 



exp 

oo 

x /« 





1 



e 4 



|<7 (s)T^(6) 



IV 



^ 4™ (1 + 717)71! 



(59) 



where 7 is defined in (25). 

As far as the Glauber exponentials are averaged over 
the eigenstates of interaction, for four different models, 
Eqs. (56)- (59), we are in a position to calculate the total, 
elastic and inelastic pA cross sections given by Eqs. (53), 
(54) and (55) respectively. 

The numerical results for four nuclei and energy y/s = 
41.6 GeV of the experiment HERA-B are displayed in 
Table I in parentheses for four different combinations of 



the models for the valence quark wave function of the 
proton and the dipole cross section. 



TABLE I: Cross sections of coherent reactions calculated 
within the Glauber model and also corrected for inelastic 
shadowing which either excludes (in parentheses), or involves 
the corrections for gluon shadowing. 



Obs. 


Nucl. 


Glauber 
model 


Model 
I 


Model 
II 


Model 
III 


Model 
IV 




W 


3073 


(2462) 
2382 


(2727) 
2641 


(2908) 
2831 


(3000) 
2918 


pA 
w tot 

[mb] 


Ti 


1159 


(938) 
915 


(1028) 
1003 


(1103) 
1077 


(1130) 
1102 




Al 


726 


(594) 
582 


(647) 
632 


(692) 
677 


(707) 
690 




C 


380 


(349) 
344 


(372) 
366 


(363) 
357 


(369) 
361 




W 


1196 


(748) 
695 


(931) 
867 


(1061) 
1001 


(1137) 
1071 


pA 

i 

" el 

[mb] 


Ti 


352 


(217) 
205 


(268) 
252 


(313) 
296 


(331) 
312 




Al 


196 


(122) 
116 


(149) 
141 


(174) 
165 


(183) 
173 




C 


79.1 


( 51.6) 
49.6 


(61.5) 
58.9 


(70.5) 
67.4 


(73.3) 
70.0 




W 




(153) 
156 


(82.3) 
86.5 


(57.2) 
59.0 


(19.2) 
20.6 


pA 

j 

sd 

(PPM) 


Ti 




(66.3) 
65.7 


(39.4) 
39.5 


(23.5) 
23.1 


(9.4) 
9.4 


[ mbl 

L J 


Al 




(42.4) 
41.8 


(25.7) 
25.5 


(14.4) 
14.0 


(6.0) 
6.0 




C 




(22.8) 
22.4 


(14.0) 
13.8 


(6.9) 
6.7 


(3.1) 
3.0 




W 


1.2 


(4.1) 
3.7 


(2.8) 
2.7 


(2.1) 
2.0 


(1.6) 
1.6 


pA 
°sd 

(3P) 


Ti 


0.6 


(1.0) 
0.9 


(0.8) 
0.8 


(0.7) 
0.6 


(0.6) 
0.6 


[mb] 


Al 


0.3 


(0.5) 
0.5 


(0.4) 
0.4 


(0.4) 
0.4 


(0.4) 
0.4 




C 


0.1 


(0.2) 
0.2 


(0.2) 
0.2 


(0.2) 
0.2 


(0.2) 
0.2 



Our predictions for p — Au collisions at RHIC and for 
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p — Pb collisions at LHC are depicted in Tables II and III 
respectively. 

As one could have expected, the nonrealistic models 
I and II, which grossly overestimate single diffraction, 
also overestimate the inelastic corrections. On the other 
hand, the more realistic models III and IV lead to quite 
moderate deviation from the Glauber model expecta- 
tions. However, this is not the end of the story, since 
the inelastic shadowing related gluonic excitations is still 
to be included. 

Note that the result for single diffraction in Tables I, 
II, III demonstrate much higher sensitivity to inelastic 
corrections than for elastic scattering. This is easy to 
understand, in the black disc limit the elastic cross sec- 
tion reaches its maximum, ttR\, and cannot be varied by 
any corrections. At the same time, diffraction vanishes, 
therefore is extremely sensitive to the transparency of the 
nucleus, it is maximal for model I, but minimal for IV. 



Diffractive excitation of the gluonic degrees of 
freedom 



The higher Fock components of the proton wave func- 
tion contain more partons, gluons, sea quarks. Appar- 
ently, the more degrees of freedom, the larger are the 
Gribov corrections. In particular, diffractive gluon radi- 
ation, corresponding to the triple-Pomeron term, makes 
nuclear matter more transparent, i.e. leads to a further 
reduction of the total pA cross section. These corrections 
will be taken into account in Sect. IX. 

Diffractive gluon radiation also contributes to the 
single diffractive process pA — > XA and creates the 
triple-Pomeron tail 1/M 2 in the effective mass distribu- 



tion. Correspondingly, the single-diffraction cross section 
Eq. (55) must be corrected for this excitation channel. 
The cross section of coherent gluon radiation on a nucleus 
is given by a straightforward generalization of Eq. (28), 



da(pA -» XA) 



dxf 



3P 



<? h / e -i* 8l (r (> a)Tj}(6)y 

4tt(1 - x F ) I \ / 



/ 



x / d 2 r T 



1 _ e -h°(r T ,s)T A (b) 



,(60) 



with the same notations. Here the first factor implies 
the absorptive corrections, analogous to those in (31). 
They are, however, much stronger than in pp collisions 
and almost terminate diffraction on heavy nuclei, except 
at the very periphery This factor is averaged over the 
proton size weighted with its wave function squared, as 
in (14). Therefore it depends on the models for the pro- 
ton wave function considered above. Since we employ the 
phcnomcnological cross section fitted to data, the possi- 
bility of gluon radiation during propagation of the proton 
though the nucleus is included. 

The integral over tt in (60) takes into account the 
diffractive channels containing gluons coherently radi- 
ated by a valence quark propagating through a nucleus 
of thickness TA{b). We make use of the smallness of the 
mean quark-gluon separation and add up incoherently 
the gluons radiated by different projectile valence quarks. 
Besides, the two models for the dipole cross section under 
discussion are almost identical at such small separations, 
r <~ 0.3 fm. 



TABLE II: Predictions for RHIC for different LRG channels. All the cross sections are for proton-gold collisions and are in 
mb. The cross sections shown in parentheses are corrected for inelastic shadowing related only to valence quark fluctuations, 
while the bottom numbers are also corrected for gluon shadowing. The cross sections which turned out to be smaller than the 
numerical accuracy of the calculations are put equal to zero. 



Model 


pPb 


pPb 


pPb] 
a sd 

L so J PPM, 


pPb] 

° sd \ 3 p 


pPb 


pPb 

cr j 

qsd 


pPb 
a tsd 


pPb 
°dd 


Glauber 


3616.8 


1446.8 




5.1 


98.6 




42.3 




III 


(3524.4) 
3457.5 


(1367.6) 
1313.8 


(31.8) 
33.3 


(7.3) 
7.6 


(95.8) 
96.2 


(3.1) 
3.1 


(41.2) 
41.4 


(3-1) 
3.1 


IV 


(3582.0) 
3514.1 


(1417.7) 
1362.3 


(8.1) 
8.9 


(5.8) 
6.3 


(98.4) 
98.9 


(0.6) 
0.6 


(42.3) 
42.53 


(0.6) 
0.6 



r 



Therefore for the sake of simplicity we use the simple proximation is greatly enhanced by the cut off imposed 
quadratic form, Uq q cx r\. The accuracy of this ap- on large separations by the first exponential factor in 
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TABLE III: The same as in Table II, but for proton-lead collisions at LHC. 



Model 


pPb 
°tot 


pPb 
el 


pPb] 
L SO IP PR 


pPbl 


pPb 


pPb 
qsd 


pPb 
a tsd 


pPb 
°dd 


Glauber 


4241.5 


1794.9 




28.8 


141.43 




22.9 




III 


(4222.9) 
4194.2 


(1778.7) 
1755.6 


(5.5) 
5.8 


31.8 
33.4 


(142.1) 
141.8 


(0.0) 
0.0 


(23.0) 
23.0 


(0.0) 
0.0 


IV 


(4235.2) 
4207.1 


(1789.9) 
1767.3 


(0.8) 
0.9 


29.5 
31.2 


(142.8) 
142.5 


(0.0) 
0.0 


(23.1) 
23.1 


(0.0) 
0.0 



Eq. (60), which gets the form, 



da(pA -> XA) 



a., 



dXF 

4<T3 9 (n)T»^ 



3iP 
2 



/ 



In 



ir 2 (l-x F ) 

, , 4(M) 



d 2 b 



l + 2e A {b, s) 



where 



e A (b, s) 



9r 2 



16i? 2 ( S ) 



K(s)a (s)T A (b) 



(61) 



(62) 



The absorptive corrections, given in (61) by the expo- 
nential factor, are factorized from the coherent diffrac- 
tive gluon radiation, given by the logarithmic factor. 
For heavy nuclei these absorptive corrections are much 
stronger than for pp collisions (see Sect. HID). They 
practically eliminate diffraction on nuclear targets except 
the very periphery. Therefore, the cross section Eq. (62) 
varies as 

Since the main contribution comes from peripheral col- 
lisions where the projectile proton finds very few nuclc- 
ons, the absorptive corrections for NN scattering consid- 
ered earlier in Sect. HID may be important and have to 
be included into (62). This is done in the same way as 
in Sect. HID, namely using the relation Eq. (31) which 
results in the suppression factor K{s) in Eq. (62), 



K ( s ) = l -TZ 



1 



4n B%{s) 



(s) 



b pp 

^ZP 



Note that although this factor decreases with energy, it 
always remains positive. Indeed, K > even in the 
Froissart limit, where [42] a to t = kirap e In 2 (s/M, 2 ) and 
B el = a' P e In 2 (s/M 2 ), where e = a P (0) - 1 w 0.08. 
The triple Pomeron slope B^p depends on xp according 
to (29). We evaluate it at the mean value of Xf- The 
magnitude of B^ p (xf) varies very slowly with energy, as 
a double logarithm. It rises only by 10% between the 
energies of HERA-B and LHC. 

Only the absorptive correction factor in (62) depends 
on the model for the proton wave function and on the 
shape of the dipole cross section. It can be evaluated 
either within the Glauber approximation, or including 



the inelastic shadowing corrections calculated in accor- 
dance with either of the four models considered above. 
We fix a s = 0.4 as was explained in Sect. HID and in- 
tegrate over xf from 0.9 up to 1 — Mq/s. The results 
are depicted in Table I for the energy of HERA-B, and 
in Tables II and III for the energies of RHIC and LHC 
respectively. 

The coherent triple-Pomeron diffraction on nuclei at 
the energy of HERA-B turns out to be amazingly small. 
This can be understood as follows. At this energy the 
value of e(s) in (30) is sufficiently small to expand the 
logarithm up to the first order of e 2 . The same is true for 
Eq. (62) provided that the absorptive corrections elimi- 
nate the contribution of large T A {b). Then, comparing 
the cross section Eq. (30) integrated over pr, and the one 
presented in Eq. (62), we see that they are related via the 
replacement, 



ppp 

-°3iP 



A(T A ) , 



(64) 



where the mean nuclear thickness for this process is given 
by 



(T A ) = jjd 2 b (e-^(n)T»^ 2 T 2 (fe) _ 



(65) 



(63) This can be estimated as, 



(Ta) 



2irwR A 



(66) 



where R A w 1.12 fm x A 1 / 3 and w ~ 0.5 fm are the nu- 
clear radius and edge thickness for the Woods-Saxon nu- 
clear density [30]. Thus, the replacement Eq. (64) leads 
to a reduction of the cross section by about a factor 
of five for heavy nuclei, A ~ 200, and much more for 
light nuclei. This explains the smallncss of the coherent 
diffractive gluon radiation by extremely strong absorp- 
tive corrections. In the black disk limit with vanishing 
edge thickness (w — ► 0) no diffractive gluon radiation is 
possible. 
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VII. QUASIELASTIC 
SCATTERING/EXCITATION OF THE 
PROJECTILE, pA -> p(X)A* 



The simplest channel of nuclear excitation is a 
quasielastic breakup of the nucleus, caused by elastic 
scattering, pN — > pN, of the projectile proton on one 
of the bound nucleons. In this case the nucleus breaks 
up into fragments without particle production. It is diffi- 
cult to control this condition experimentally, but is easy 
to calculate it. 

One should modify Eq. (46), first sum the nuclear ex- 
citations A* and integrate over impact parameter s, then 
average over the quark positions r, and in each of the 
two amplitudes, 



a^(pA^pA*) 



= J d2& |E(°l rP ' 4 ( 6 '^)|^)I,n( F l rP ' 4 ( 6 '^)| )^ 

- (o|r^(6,r;)|o); iri (o|r^(6,^)|o) 5ir , 
= J d 2 6|(o|[r" A (6,F i )] + r^(6,fj.) 

- <o|r p >,ri)|o)* ri } . 



(67) 



Here s is the impact parameter of the projectile proton 
relative to the bound nucleons, and fi are the relative 
transverse positions of the valence quarks in the projectile 
proton. After integration over impact parameter s [see 
in (46)] we get, 



x < exp 



k=l 



k\ 



47rS e( 



d{T A f 



exp 



(68) 



This series quickly converges due to smallness of the elas- 
tic cross section. Even the fist term provides a reasonable 
accuracy. We control the accuracy to be within 1%. 

If we sum over all excitations of the proton contain- 
ing no radiated gluons and apply the condition of com- 
pleteness, a delta- function 5(fi — r[) emerges, leading to 
the following expression for the sum of quasielastic and 



quasi-single-diffractive channels , 

a^(pA^pA*) + a p q t d (pA 

d 2 b(cxp [-a(fi)T^(b)] 



XA*) 



exp 



167TB, 



el 



"Eh 

fe=i 



()'■ 



2 k 



d(T A , 



21: 



(exp[-a(rl)T>)]) n 



(69) 



With these expressions, Eqs. (68)-(69), we can calcu- 
late the quasielastic and quasidiffractive cross sections 
employing different models for the dipole cross section 
and the proton wave function. 
Model I 

In this simplest case of the mesonic type wave function 
of the proton and a(fi) = <Jq q (rT) oc r\, the quasielastic 
and quasi-diffractive cross sections read respectively, 

k 



pN 



T h A {b) 



1 + I^T» 



pA 
J qel 



pA 



E 

k=l 



(2fc)! 
k\ 



d 2 b- 



pN 



2k+2 ■ 



n(b) 



(70) 



2fe+l ■ 



(71) 
Model II 

In the case of a symmetric valence quark wave function 
and dipole cross section a(fi) = as q (fi) the quasielastic 
cross section takes the form, 



pA 
'qel 



E 

fe=l 



(fc + l)(fc + 1)! 

2 2k 



pN 



n(b) 



1 + %<x. 



2k+4 ' 



>t mb) 

(72) 

This result may look surprising. Indeed, the quasielas- 
tic cross section Eq. (47) is proportional to nuclear trans- 
parency, i.e. the survival probability of a proton propa- 
gating through the nucleus. That is given by the mean 
value Eq. (57) squared, i.e. the fourth power of T A in 
the denominator. That would mean more nuclear trans- 
parency and larger quasielastic cross section compared to 
the Glauber model. Eq. (72), however, has the leading 
term which behaves as the inverse sixth power of T A . Al- 
though at large T A the exponential Glauber transparency 
is less than any power of T A , it turns out that for real 
nuclei inelastic shadowing makes nuclei less transparent, 
at variance with the simplified expectation. The reason 
why the nucleus is less transparent than is suggested by 
Eq. (57) is an extra weight factor <Jq q {r) in the quasielas- 
tic amplitude. This factor suppresses small-size projectile 
components for which the nucleus is transparent. 
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Using completeness we can also calculate the sum of 
the cross sections of quasielastic and quasidiffractive scat- 
tering, the result reads, 



pA . pA 
^qel ' ®qsd 



£ 

k=i 



(2fe+l)l 
2 2k k\ 



d 2 b- 



pN 



TH(b) 



2k+2 ' 



(73) 

In this case the leading term behaves like T A 6 at large 
Ta, i.e. very heavy nuclei are much more transparent for 
quasidiffractive, than quasielastic processes. 
Models III and IV 

We skip the cumbersome expressions for the quasielastic 
and quasi-diffractive cross sections in this case. We use 
instead equations (68) and (69) respectively and perform 
calculations numerically. For the averaged Glauber ex- 
ponential (exp [— \(r{fi) T^(6)] ) we employ Eqs. (58) 
and (59) for models III and IV respectively. 

So far we calculated the quasi-diffractive cross section 
related to proton excitations without gluon radiation. 
Now we should include diffractive gluon brcmsstrahlung, 
as it is done for coherent diffraction in Sect. VI B. This 
can be done via a simple replacement in the above equa- 
tions (68)-(73), 



pA 



a 



The results for cr^ and o v * d are depicted in paren- 
theses in Table IV for different nuclei at the energy of 
HERA-B, and in Tables II and III for gold at RHIC and 
lead at LHC respectively. In both models III and IV the 
quasi-diffractive and double-diffractive cross sections are 
very small due to very low nuclear transparency close to 
the black disc limit, so small at the energy of LHC that 
in model IV we could not reach a sufficient numerical ac- 
curacy. Therefore these cross sections are set to zero in 
Table III. 

One can notice a much higher sensitivity to the models 
for quasidiffractive scattering compared to quasielastic. 
In the former case this simply reflects the tremendous 
model dependence of single diffraction, pointed out in 
Sect. III. In the latter case there is a partial compensa- 
tion between the steep rise of the dipole cross section with 
the dipole size (oc r\ at small Tt) and the exponential 
fall of nuclear transparency. 



sdipPM 
pA 



pA 
qsd 



(74) 



VIII. DIFFRACTIVE EXCITATION OF BOUND 
NUCLEONS 

Another LRG process, not explored so far, is diffrac- 
tive excitations of bound nucleons in the target. It may 
happen along with or without excitation of the beam. 
We call these cases target single diffraction, or double 
diffraction respectively. 

Since the phenomenological dipole cross is averaged 
over the target proton wave function, it is insufficient 



for these reactions which need knowledge of a doublc- 
dipole cross section. If, however, one assumes that in pp 
collisions the dependence on the beam dipole size is the 
same for reactions with or without target excitation, then 



TABLE IV: Same as in table I, but involving diffractive break 
up of the nucleus. 



Obs. 


Nucl. 


Glauber 
model 


Model 
I 


Model 
II 


Model 
III 


Model 
IV 




W 


88.2 


(59.1) 
59.0 


(45.9) 
46.0 


(73.0) 
73.5 


(77.3) 
78.0 




Ti 


63.9 


(42.5) 
42.0 


(35.3) 
35.2 


(53.2) 
52.8 


(56.0) 
55.5 




Al 


48.8 


(32.8) 
32.4 


(30.0) 
27.8 


(40.7) 
40.3 


(42.7) 
42.2 




C 


34.5 


(23.2) 
22.8 


(20.5) 
20.3 


(28.6) 
28.1 


(29.7) 
29.2 




W 


- 


(41.9) 
40.7 


(22.6) 
22.3 


(6.2) 
6.1 


(2.5) 
2.5 




Ti 


- 


(31.4) 
30.6 


(18.7) 
18.6 


(4.4) 
4.5 


(1.8) 
1.9 




Al 


- 


(24.8) 
24.2 


(15.5) 
15.4 


(3.5) 
3.6 


(1.5) 
1.5 




C 


- 


(18.2) 
17.7 


(11.6) 
11.6 


(2.7) 
2.7 


(1.2) 
1.2 




W 


42.1 


(28.2) 
28.1 


(21.9) 
22.1 


(34.8) 
35.1 


(36.9) 
37.2 


Ob] 


Ti 


30.5 


(20.3) 
20.1 


(16.9) 
16.8 


(25.4) 
25.2 


(26.7) 
26.5 




Al 


23.3 


(15.7) 
15.5 


(13.3) 
13.3 


(19.4) 
19.2 


(20.4) 
20.1 




C 


16.5 


(11.0) 
10.9 


(9.8) 
9.7 


(13.6) 
13.4 


(14.2) 
14.0 




W 




(15.0) 
14.5 


(8.1) 
8.0 


(2.2) 
2.2 


(0.9) 
0.9 


rfd [ mb ] 


Ti 




(11.2) 
10.8 


(6.7) 
6.7 


(1.6) 
1.6 


(0.7) 
0.7 




Al 




(8.9) 
8.6 


(5.5) 
5.5 


(1.3) 
1.3 


(0.5) 
0.5 




C 




(6.5) 
6.3 


(4.3) 
4.1 


(1.0) 
1.0 


(0.4) 
0.4 
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one can relate the cross sections in question to known 
quasielastic and quasidiffractive ones, 

pp 

^PY) = -§£ a^ipA pA*) ; (75) 

"el 
-PP 

o£{pA^XY) = JL^ d { P A^XA*): (76) 

Numerical results are depicted in Table IV for HERA-B 
and in Tables II and III for RHIC and LHC respectively. 

IX. GLUON SHADOWING AND THE 
TRIPLE-POMERON DIFFRACTION 

Eikonalization of the lowest Fock state \iq) of the pro- 
ton done in (53) corresponds to the Bethe-Heitler regime 
of gluon radiation. Indeed, gluon bremsstrahlung is re- 
sponsible for the rising energy dependence of the phe- 
nomenological cross section (3), and in the eikonal form 
(36) one assumes that the whole spectrum of gluons is 
multiply radiated. However, the Landau-Pomeranchuk- 
Migdal (LPM) effect [43, 44] is known to suppress radia- 
tion in multiple interactions. Since the main part of the 
inelastic cross section at high energies is related to gluon 
radiation, the LPM effect becomes a suppression of the 
cross section. This is a quantum-mechanical interference 
phenomenon and it is a part of the suppression called 
Gribov's inelastic shadowing. The way how it is taken 
into account in the QCD dipole picture is by the inclu- 
sion of higher Fock states, \3qG), etc. Each of these states 
represents a colorless dipole and its elastic amplitude on 
a nucleon is subject to eikonalization. 

As we already mentioned, the eikonalization procedure 
requires the fluctuation lifetime to be much longer than 
the nuclear size. Otherwise, one has to take into ac- 
count the " breathing" of the fluctuation during propaga- 
tion through a nucleus, which can be done by applying 
the light-cone Green function technique [1, 34, 45]. In 
hadronic representation this is equivalent to saying that 
all the longitudinal momenta transfers must be much 
smaller than the mean free path of the hadron in the 
nucleus. Otherwise, the phase shifts and interferences 
are to be included [46] . 

The cm. energies of HERA-B, RHIC and LHC are 
sufficiently high to treat the lowest Fock state containing 
only valence quarks as "frozen" by the Lorentz time di- 
lation during propagation through the nucleus. Indeed, 
for the excitations with the typical nucleon resonance 
masses, the coherence length is sufficiently long compared 
to the nuclear size. This is why we applied eikonalization 
without hesitation so far. Such an approximation, how- 
ever, never works for the higher Fock states containing 
gluons. Indeed, since the gluon is a vector, the integra- 
tion over effective mass of the fluctuation is divergent, 
dM 2 /M 2 , which is the standard triple-Pomeron behav- 
ior. Therefore, the colliding energy never can be suffi- 
ciently high to neglect the large-mass tail. For this rea- 



son the inelastic shadowing corrections related to exci- 
tation gluonic degrees of freedom and calculated in the 
tree-approximation never saturate, and keep rising loga- 
rithmically with energy. 

There are, however non-linear effects which are ex- 
pected to stop the rise of inelastic corrections at high 
energies. This is related to the phenomenon of gluon 
saturation [47, 48] or color glass condensate [49]. The 
strength of these nonlinear effect is model dependent. It 
is a rather mild effect within the model of small gluonic 
spots explained in Sect. Ill C. The reason is simple, in 
spite of a sufficient longitudinal overlap of gluon clouds 
originated from different nuclcons, there is insufficient 
overlap in the transverse plane. This fact leads to a de- 
lay of the onset of saturation up to very high energies, 
since the transverse radius squared of the gluonic clouds 
rise with energy very slowly, logarithmically, with a small 
coefficient of the order of 0.1 GeV~ 2 

The details of the calculation of inelastic corrections 
related to excitation of gluonic degrees of freedom can 
be found in [1, 16]. The results for nuclear cross sections 
corrected for gluon shadowing are shown in Tables I, II, 
III, and IV. 



X. CONCLUSIONS 

In this paper we study the dynamics of LRG processes 
in pp and pA collisions basing on the light-cone dipole ap- 
proach which allows to develop phenomenology not only 
for hard, but also for soft reactions. 

One of our objectives is an improvement of accuracy of 
calculations for nuclear effects in LRG processes. First of 
all, we tested some popular models for the dipole-proton 
cross section and the proton wave function on their abil- 
ity to reproduce the cross section of diffractive excitation 
of the valence quark skeleton of the proton. Data show 
that the forward diffractive cross section is amazingly 
small, only about 10% of the elastic one. We conclude 
that the models based on the dipole cross section which 
quadratically depends on quark separation fail badly to 
explain even the order of magnitude of the single diffrac- 
tive pp cross section. This model, however, is quite pop- 
ular in the literature devoted to nuclear effects, in par- 
ticular color transparency, since it helps to simplify the 
calculations. Apparently, those results cannot be realis- 
tic. On the other hand, we found the saturated shape of 
the dipole cross section which levels off at large separa- 
tions to be quite successful in explaining the data on pp 
diffraction. 

As for nuclear effects, the popular Glauber model can- 
not treat properly most of the off-diagonal diffractive 
processes, since this is a single-channel approximation. 
Based on the color-dipole representation we develop tech- 
niques for calculating the cross sections of LRG processes, 
both diagonal and off-diagonal. This method allows to 
sum the Gribov inelastic shadowing corrections to all or- 
ders. These corrections make nuclear matter more trans- 
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parent and reduce the total hadron-nucleus cross section. 
At the same time, their influence on other diffractive 
reactions depends on a complicated color transparency 
interplay making nuclei more transparent for small size 
hadronic fluctuations, but simultaneously suppressing 
the strength of the interaction with bound nucleons, due 
to the same effect. This is confirmed by our numerical 
results for the cross sections of a variety of channels pre- 
sented in Tables I-IV. We found that models I and II 
based on the quadratic dependence of the cross section 
on the dipole size grossly overestimate the Gribov cor- 
rections compared to more realistic variants III and IV, 
based on the saturated form of the dipole cross section. 

Available data for the Gribov corrections [40, 41] show 
that they rise with energy, what results from the falling 
energy dependence of the longitudinal momentum trans- 
fer. At sufficiently high energies the PIPM part of the 
single diffraction saturates, since the longitudinal mo- 
mentum transfer vanishes, while the triple-Pomeron part 
keeps rising logarithmically. Comparing, however, the re- 
sults depicted in Tables I and III we see that the inelastic 
corrections, i.e. the deviation from the Glauber model, 
vary from 5 — 7% at the energy of HERA-B down to about 
1 — 2% at LHC. Such a dramatic reduction signals the 
closeness of the unitarity limit. Indeed, in this regime, 



also called black-disk limit, different Fock components of 
the proton interact with the same cross section. There- 
fore the incoming and outgoing waves consist of the same 
superposition of Fock states, then only elastic scattering 
is possible. Similar suppression of diffraction is expected 
for pp scattering when it reaches the Froissart regime, 
a sd/ a tot ^ l/hi(s). However, the onset of this effect on 
nuclear targets can be observed at much lower energies. 
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